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We study quantitatively the formation and evolution of be, be bound states in a space-time domain 
of deconfined quarks and gluons (quark-gluon plasma, QGP). At the Relativistic Heavy Ion Collider 
(RHIC), one expects for the first time that typical central collisions will result in multiple pairs of 
heavy (in this case charmed) quarks. This provides a new mechanism for the formation of heavy 
quarkonia which depends on the properties of the deconfined region. We find typical enhancements 
of about 500-fold for the fee, 6c production yields over expectations from the elementary coherent 
hadronic _Bc-meson production scenario. The final population of bound states may serve as a probe 
of the plasma phase parameters. 
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I. INTRODUCTION 

The recent observation of candidate iJc-meson events 
by the CDF collaboration yields measurements for 
the ground state mass and lifetime which are consistent 
with expectations from nonrelativistic potential models 
1^ . This system of states is expected to have properties 
intermediate between the J /ip and T systems (except for 
its longer lifetime due to the absence of quark annihi- 
lation processes). Thus it is natural to investigate the 
fate of such states when produced in a deconfined envi- 
ronment, where suppression relative to "expected" yields 
may serve as a signature for the existence of deconfine- 
ment. 

We have previously reported some preliminary work 
toward this end |^,^. Our calculations of the expected 
yield of Be at RHIC utilize estimates of the coherent 
production of both a bb and cc pair in the same initial 
hard perturbative QCD process, followed by hadroniza- 
tion into the Be states. This process is of order a^, and 
falls quite rapidly with decreasing energy. The result- 
ing prediction at RHIC energies for the bound-state frac- 
tion relative to initial b-quark production is in the range 
10~^ — 10~^ [^. Combining this result with the expected 
yield of b-quarks, we find that even at design luminosity 
there will be at most a handful of Sc-mesons produced at 
RHIC which decay into observable final states (all scalar, 
vector, excited, particle/antiparticle be bound states will 
be called i3c-mesons or simply B^ except when otherwise 
noted). Thus a scenario of production via calculable hard 
QCD interactions followed by suppression in a deconfined 
medium will be irrelevant for RHIC parameters. 

What we now present is a new production mechanism 
which itself depends on the existence of a deconfined 
state. The new mechanism becomes relevant for the 
first time at RHIC energies, since typical central colli- 
sion events will have multiple pairs of initially-produced 
charm quarks. Estimates using perturbative QCD to cal- 
culate the initial production of heavy quarks predict ap- 
proximately 10 cc pairs in each central event in Au-Au 



collisions at ^/s = 200 GeV Due to the larger mass, 
only about 1 in 20 central events will produce a bb pair. 

Consider the events in which there is a single bb pair, 
along with the expected 10 cc pairs. If a region of decon- 
fined matter is subsequently produced in the space-time 
region encompassing the heavy quarks, a b-quark will be 
able to "find" any of the initially-produced charm quarks 
with which to produce the final Be bound state. More 
generally, our study shows that the formation of quarko- 
nium states within a deconfined space-time domain re- 
sulting from the interaction of mobile heavy quarks offers 
a very interesting signature for deconfinement. 

In Sec. 1^ we study quantitatively the probability that 
in the deconfined quark-gluon plasma (QGP) phase the 
mobile c, c-quarks can seek, find and bind to a 5, 5 quark 
produced in the same event. We also consider the dy- 
namical evolution of these h eavy quark bound states in 
the deconfined phase. In Sec. [I A we establish the kinetic 
model for the evolution in time of the be, be bound state 
popula tion . The required cross sections are presented 
in Sec. 



II B 



and the formation and dissociation reaction 
rates obtained. 

To calculate the final Be yield, one needs an estimate 
of the charm quark density during the deconfined phase. 
We therefore consider in Sec. [II A a generic model for the 
expansion and cooling of the QGP. This also enables the 
comparison of direct (i nitial) and microscopic thermal 



charm production in Sec. [II B. We confirm that at RHIC 



energy charm is produced primarily in the initial parton 
interactions, with only a minor contribution arising from 
the thermal plasma processes. We also show that the 
charm density at hadronization remains significantly in 
excess of that which would be present if full chemical 
equilibrium were reached. 



We then present in Sec. [V a study of the influence of 



the model parameters on the pattern of Be production 
at RHIC. We present in Sec. IVB the fractional yield per 



initial bottom quark pair in QGP. We find that the fi- 
nal yield is most sensitive to the initial charm density, 
and there is very little dependence on the initial temper- 



1 



ature o f the dense deconfined state. In contrast, we show breakup reaction rate has its most significant effect [O 



in Sec.|VD| that the relative yield of the first radially 
excited state to the ground state, Bc{2S)/Bc{lS) is sig- 
nificantly more dependent on the initial temperature in 
the plasma phase. 



A. Chemical rate equations 



II. CHEMICAL KINETICS OF QUARKONIUM 
ABUNDANCE 

In our scenario, the be and be states formed in the 
QGP will be subject to collisions with gluons, which 
will dissociate them into their constituent quarks. This 
mechanism can be thought of as the dynamic coun- 
terpart of the plasma screening scenario, in which the 
color-confinement force is screened away in the hot dense 
plasma |||,^ . We do not distinguish between the vacuum 
and plasma values for the mass and binding energy of the 
IS ground states. Both are significantly larger than typ- 
ical temperatures expected at RHIC. Estimates for their 
behavior as a function of screening mass have been made 
and indicate that such an approximation is reason- 
able. Thus we no longer distinguish between the plasma 
bound states and the physical Be mesons which can be 
observed after hadronization. 

The primary formation mechanism is just the inverse 
of the breakup reaction, in which unbound heavy quarks 
are captured in the bound states, emitting a color octet 
gluon. The competition between the rates of these reac- 
tions integrated over the lifetime of the QGP then deter- 
mines the final Be populations. Note that in this scenario 
it is impossible to separate the formation process from 
the breakup (suppression) process. Both processes oc- 
cur simultaneously, in contrast to the situation in which 
the formation only occurs at the initial times before the 
QGP is present. (Of course, one can include this case 
also in our scenario by adjusting the initial conditions. 
However, for the case of Be initial production at RHIC, 
we have already shown that this possibility is negligible.) 
A number of other reactions involving b-quarks are pos- 
sible in the QGP, but the rates are much smaller than 
those above. For example, formation of bound states 
with light quarks, Bg, Bu, Bd, are not possible at the high 
initial temperatures expected at RHIC since their lower 
binding energies prevents them from existing in a hot 
QGP, or equivalently they are ionized on very short time 
scales. They will be formed predominantly at hadroniza- 
tion, but this process is too late to affect the final Be pop- 
ulation. We also neglect the decrease in b-quark popula- 
tion due to bb annihilation into light quarks or formation 
of T states. Both of these processes depend quadratically 
on the b-quark densities, and proceed too slowly or too 
rarely, respectively, to be significant. We also neglect the 
breakup cross section due to collisions with light quarks, 
since their population in the QGP is expected to be very 
much suppressed relative to gluons as compared to chem- 
ical equilibrium values during the early times when the 



The abundance of bottom and charm quarks and 
their bound states is thus governed by very simple mas- 
ter (population) equations involving only two reactions. 
Specifically, the formation reaction F 

b + c^ Be + g , 
and the dissociation reaction D 



Be 
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b + c, 



and similar equations for the conjugate states determine 
the time evolution of the number of bound states in the 
deconfined region, Nb^ ■ 



dNB. 
dT 



= XpNb Pc - XbNb, Pg . 



(1) 



Since the total number of b-quarks does not change under 
the assumptions above, the rate of change in the number 
of unbound b quarks Nh is just the negative of that for 
the number of bound state Be mesons : 



dNb 
dT 



= XdNb, Pg - Af A^6 Pc 



(2) 



In Eqs. (0j^), T is the proper time in a small volume cell. 
Pi denotes the number density [L~^] of species i and the 
reactivity A [L'^/time] is the momentum distribution av- 
eraged reaction rate: 



A = (crUrcl) = 



/ / d^pid^p2flipi)f2iP2)cr{^/l)Vr:cl 
J J d^p,d^P2fl{pi)f2{p2) 



(3) 



Note that even though A and p are not Lorentz invariant, 
their product, the rate of particle production, is invari- 
ant. Thus, we evaluate the rates PiXiy/p in the volume 
cell frame, where the (local) densities are given by exter- 
nal inputs. 

We study here the deconfined period during which par- 
ton distributions have kinetically (but not necessarily 
chemically) equilibrated. We thus consider an isotropic 
medium of mobile quarks and gluons with the momentum 
distribution function fi for particle i given by the thermal 
equilibrium distribution / — (e^/-^ ± 1)~^. cr is the spin- 
and color-averaged total cross section for the relevant re- 
action which depends only on s = {pi +P2)'^, and v^d 
is the relative speed between the two reacting particles. 
Except where otherwise indicated pg is assumed to be in 
chemical equilibrium, pe will be determined by its own 
kinetic equation described in Sec. IIIB. The spatial dis- 



tribution of charm and bottom quarks will be taken to be 
uniform, an approximation which indicates the accuracy 
to which we will pursue our following calculations. 
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B. In-plasma cross section and reaction rates 

We now need to estimate the cross sections involved 
in the formation and breakup reactions. For these pur- 
poses, we first utilize a derivation based on the inter- 
action of a gluon field with the color dipole moment of 
a nonrelativistic heavy quarkonium state. It is imple- 
mented via an operator product expansion technique , 
and has been applied to the J/il) breakup rates in a QGP 
We generalize this result to any heavy quarkonium 
IS state with arbitrary flavor content, so that the spin 
and color averaged dissociation cross section is written 



27r /32\2/2^\i/2fc^/2 



ko) 



3/2 



(4) 



Here k is the gluon momentum (in the quarkonium-rest 
frame), fcg is the minimum value required to impart the 
binding energy e to the bound state, and /i is the reduced 
mass. This form is valid if the quarkonium system has 
a spatial size small compared with the inverse of Aqcd, 
and its bound state spectrum is close to that in a nonrel- 
ativistic Coulomb potential with e large compared with 
^QCD- These conditions are marginally satisfied for the 
J/ip, and should be somewhat better for the Be kinemat- 
ics. The form above has also been altered to account for 
recoil of the finite mass system, since in the original form 
the values of e and ko were identical. 

We use values mi, = 5.8 GeV and rric = 1.3 GeV which 
are consistent with typical potential model fits to the 
spectra. The magnitude of the cross section is controlled 
by the geometrical factor (4/i^)~^, which for these quark 
masses is consistent with the size of the bound state wave 
function in the same potential models. The rate of in- 
crease just above threshold is due to phase space and the 
p-wave color dipole interaction, and it reaches a max- 
imum value for dissociation of about 1.5 mb when 
k = -y fco . For our model calculations we use a central 
value of 6.3 GeV for the Be mass, and a binding energy 
of 0.84 GeV which follows from the hadronic open fla- 
vor B and D meson masses. We retain these values for 
our kinetic calculations, as an approximation for physical 
values in the QGP. One might expect that both would 
decrease somewhat in a QGP. Some potential model cal- 
culations utilizing screening indicate the magnitude of 
this effect is expected to be small for the Be at RHIC 
conditions pC| ]. 

Our dominant formation cross section trp can then be 
directly obtained from an by utilizing the detailed bal- 
ance relation. This is written in the zero-momentum 
(ZM) frame of two-body interactions as 



(5) 



where Pd/f is the 3- momentum of the initial state par- 
ticles for these reactions in their respective ZM frames, 
and gu/F is the statistical degeneracy in the two chan- 
nels. In the present case of an unpolarized QGP, go = 



((3-1- 1) • 1)(2 • 8), counting spin and color multiplicities of 
Be color-singlet and gluon, and gp = {2 ■ 3) (2 • 3), count- 
ing spin and color multiplicities of the initial state b- and 
c-quarks. In this formalism we have included both the 
pseudoscalar and vector IS i?c states with the same cross 
sections, since in the nonrelativistic approximation one 
would expect the same spatial wave function and no spin 
dependence. These two states are also included in the no- 
deconfinement scenario initial production estimates 
so that a direct comparison can be made. Using these 
formation - dissociation cross sections we calculate the 
reactivities as defined by Eq. As shown in Fig.[|, 
formation and dissociation tend to balance at very high 
temperature, where the endotherm nature of the dissoci- 
ation is negligible. At temperatures in the range 150-300 
MeV the formation reactivity dominates by about a fac- 
tor 4. 
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FIG. 1. Formation Ap (solid) and dissociation Ad (dashed) 
reactivities as function of temperature. 



C. Be dissociation time scales 

From Eq. (|l|), one sees that if the system is at a con- 
stant temperature (or equivalently the reaction rates 
are sufficiently fast), the final ratio of bound state 
to free quark populations is just given by the ratio 
(AF/Ce)/(ADpg). The relevant time scales are set by the 
magnitudes of either factor in this ratio. We have calcu- 
lated the dissociation rates of several quarkonium states 
under breakup by gluons with full chemical equilibrium 
density, and the results are shown in Fig. ^j. 

As expected, these rates rise quite sharply with tem- 
perature. The Be curve fits quite nicely between those 
for the J/ip and T, also as expected. We also show the 
same calculation for the Bs state, although the approxi- 
mations made for this cross section have a very marginal 
validity in view of such a large state with small binding 
energy. For the Be, one sees that in the range of ini- 
tial temperatures expected at RHIC (roughly 300 to 500 
MeV) , these dissociation rates imply time scales of order 
1-10 fm/c. Since the total QGP lifetimes are also in this 
region, it is evident that the equilibrium solutions will 
not in general be reached at each temperature, and one 
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must solve the rate equations numerically to obtain the 
final populations. 
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FIG. 2. Thermal QGP quarkonium dissociation rates by 
thermal gluons as functions of temperature. 



III. EVOLUTION OF CHARM DENSITY IN QGP 
A. Evolution of temperature 



Note that we have rescaled the initial proper time to zero. 

This simple approach produces temperature vs. time 
profiles which appear to be very similar to those aris- 
ing in more complex studies. These are shown in Fig.^ 
for a homogeneous bulk QGP state in a range of initial 
temperatures 0.3 < To < 0.6 GeV. Intersection with the 
horizontal line at freeze-out temperature T/ — 0.15 GeV 
indicates the QGP lifespan between thermalization and 
hadronization. 
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FIG. 3. Temperature [GeV] vs. proper time [fm] for var- 
ious initial temperatures To, with parameters: li = \ fm, Vr 
— 0.58. The horizontal line marks hadronization. 



In order to obtain the chemical evolution of quark 
bound state abundances in QGP we need to establish 
a relation between plasma temperature and proper time. 
We assume that the expansion of the QGP follows an 
isentropic path We utilize a generic scenario for 

the proper-time dependence of the volume involving both 
longitudinal and transverse expansion, and examine the 
sensitivity of our final results to variations in the param- 
eters involved. As a specific example we consider an adi- 
abatically expanding homogeneous QGP domain having 
the shape of an ellipsoid of revolution about the longitu- 
dinal axis with semi-major and semi-minor axes param- 
eterized with an initial length li/2 = tq and an initial 
transverse radius r^. These are fixed at the time of QGP 
equilibration at an initial temperature Tq, when our for- 
mation and breakup reactions are assumed to start. We 
will explore the range 0.5 < U < 2 fm (equivalent to ther- 
malization times between 0.25 and 1.0 fm), and take = 
5 fm for Au-Au collisions at RHIC, corresponding to 15- 
20% most central interactions. The longitudinal growth 
occupies the region between the (almost unstopped) re- 
ceding nuclei. For transverse growth we allow a radial 
expansion at a speed Vr- The speed of sound in an ideal 
relativistic gas Vr — 0.58 c is taken as a nominal value, 
but the effects of significant variations in parameter space 
will be considered. Thus the volume evolves as a function 
of proper time r according to: 

V[r) = ^{n+VrTf{^^+T). (6) 



B. Thermal charm production in QGP 

From fits to experimental nucleon-nucleon reaction 
data it is estimated that at RHIC there will be an av- 
erage number Nc^ = 10 of directly produced cc pairs per 
central collision [^, and this is the standard value for 
which our calculations will be carried out. However, ad- 
ditional charmed quarks could be produced in the QGP 
by collisions of gluons: 

g + g -> c + c, 

and of light quarks: 

q + q-^ c + c. 

The local evolution of charm density pc can be described 
in a fashion quite similar to the model developed for 
strangeness production, see e.g. |l^. Allowing for charm 
production and volume dilution under adiabatic expan- 
sion, the local charm density obeys: 

where T = dT/dr . The rate of charm production is 

MT] - IpfX''^'' + pfX"'^^'' , (8) 
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where p°° denotes the density in thermal and chemical 
equilibrium. Ac[r] is shown in Fig.^, obtained with run- 
ning QCD parameters, as{Mz) = 0.118, mc(lGeV) = 
1.5 GeV |l|. 




T [GeV] 

FIG. 4. (Invariant) charm production rate per unit of time 
and volume Ac[T] as function of temperature T. 



Integrating Eq.|^, we obtain the total number of 
charmed quark pairs in plasma at freeze-out, with ini- 
tial temperature Tq being a parameter, along with the 
initial longitudinal size h . Fig. ^ shows our results for 
three different initial longitudinal sizes li = 2fm (dot- 
ted), k = Ifm (dashed), and k = 0.5fm (solid): the 
three thick up-curving lines include both the 10 directly 
produced charm pairs and the QGP charm production. 
The larger the initial volume and the longer the lifespan 
at high temperature, the greater is the QGP contribution 
to charm yield. 
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FIG. 5. Thick lines: Nc at freeze-out as function of initial 
temperature To [GeV]: thermally produced charm is combined 
with Nco = 10 charm pairs from direct initial production. 
Longitudinal sizes are h = 2 im (dotted), 1 fm (dashed), and 
0.5 fm (solid); other parameters same as Fig.0. Thin lines: 
chemical equilibrium abundance at freeze-out, V p'^lxf 

For the range of initial temperatures expected at 
RHIC, ones sees that there is virtually no additional 
charm produced in the QGP, and also that charm an- 



nihilation processes are too slow to significantly reduce 
the initially-produced number of charm quark pairs. The 
thin lines show for comparison the number of charm 
quarks which would follow from a chemical equilibrium 
density at the hadron freeze-out temperature Tf = 0.15 
GeV for rric = 1.5 GcV. (This value depends on the ini- 
tial temperature To through its effect on the freezeout 
volume.) The chemical equilibrium value is always sig- 
nificantly smaller than the corresponding direct -I- ther- 
mally produced abundance, i.e. direct initial charm pro- 
duction at RHIC is predicted in general to significantly 
oversaturate the statistical phase space at freeze-out. 

IV. DISCUSSION OF RESULTS 

A. Be survival in plasma 

To see what temperature range contributes to breakup, 
Fig. ^ shows the survival probability of one Be meson 
being placed in the QGP at different initial temperatures. 
At low initial temperature {T < 350 MeV) the state will 
most likely remain bound, while already at T > 500 MeV 
the plasma will most likely dissolve the Be- One can infer 
from this result that the final state population is mainly 
dominated by the formation process as the plasma cools 
between these temperatures. This of course is not in 
exact correspondence with the actual physical process, 
which involves the formation of the bound state at any 
time after QGP forms at the initial temperature, where 
the expansion has already produced a partial decrease in 
the gluon density. 




FIG. 6. Survival probability of an initial Be as function 
of initial temperature To [GeV]; other parameters same as 
Fig.^ Solid line: Nc — 104-thermal production, dashed line: 
iVc = 10. 



Another way to interpret the results shown in Fig.g 
is that the effective collision-mediated color screening is 
most effective in dissolving Be bound state when the av- 
erage gluon energy Eg ~ 3T exceeds the binding energy. 
Using this argument we can also easily see that the fla- 
vor transfer reaction bs ~\- c ~* be + s is unlikely, since the 
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bs bound states have much smaller binding energies and 
are color screened already at all temperatures above the 
hadronization temperature. 

B. Fractional Be yields 

We numerically integrate the evolution of the B^. abun- 
dance as function of the temperature T, using our ex- 
pansion model to obtain the temperature as a function 
of time. The initial conditions are Nb^ = 0, and Nb = 
1. We use Ncg = 10 and the volume expansion model 
for the charm quark density, and a full thermal gluon 
density. The final Be population per initial bb pair is 
shown in Fig.|^. The factor of 2 in the axis label indi- 
cates that on average an equal number of particle and 
antiparticle bound states will be produced. Thus the 
numbers can be directly compared with the bound state 
ratio 10^^ — 10^^ expected if no deconfinement occurs. 

0.07 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 
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FIG. 7. Be fractional abundance vs. temperature [GeV] 
for various To until freeze-out. Dashed line: Pg — ■ Other 
parameters same as Fig.^. 

The nominal initial volume and expansion parameters 
are used. We explore a range 0.2 < To < 0.6 GeV of 
initial temperature in steps of 0.1 GeV. We see that the 
Be abundance grows approximately linearly during the 
entire expansion as temperature decreases. This veri- 
fies that the formation and dissociation reactions do not 
come to equilibrium during the QGP phase, as previously 
anticipated in This linearity is not of any inherent 
physical origin, as can be seen by changing variables in 
Eq.|^ from r to T and noting that Nb, < N^. The lin- 
earity in T is seen to follow from the numerical values of 
the product T[t] V[t\ Af remaining almost constant for 
To ^ 0.6 GeV . The important conclusion we draw from 
the results shown in Fig.0 is that we can expect a rather 
To independent fractional Bc-yield, which for the main 
benchmark of our assumptions is at 5%. 

Shown by dashed lines in Fig.|^ is the scenario where 
gluon density is 1/2 its thermal equilibrium value p^{T). 
This corresponds to an effective reduction in the available 
degrees of freedom for gluons in a QGP for temperatures 



that are not yet 'asymptotically' large ||15|. One sees the 
expected effect of reduced gluon density leading to re- 
duced dissociation, and hence increased final bound state 
populations. However, the increase is substantially less 
than linear, indicating that the formation term in the 
rate equation is dominant. 

In order to better understand this result, and also to 
illustrate the dependence on initial QGP volume and 
transverse velocity, we show in Fig. ^ the freeze-out frac- 
tional i?c-yield per one bottom quark pair, as a function 
of To, while also varying the other parameters. 

0.08 I 1 1 1 1 1 1 1 
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FIG. 8. Be fractional abundance at freeze-out vs. To for 
Vr = 0.58c (solid), Vr = 0.79c (dashed), Vr = c (dotted) and 
Vo = 26 fm'^ (upper 6 lines), 52 fm"^ (middle 6 lines), 105 fm"^ 
(lower 6 lines). The effect of neglecting additional charm pro- 
duction in the QGP is shown in the thin lines which are only 
visible for T > 0.5 GeV. 

As To increases, the yield initially increases. However, 
as the initial temperature further increases, the disso- 
ciation reaction becomes more prominent and adversely 
impacts the yield, an effect which is not fully compen- 
sated by the additional QGP-produced charm available 
at these temperatures. This is seen comparing the thick 
up-curving lines with thin lines which do not include 
QGP-produced charm. However, this competition be- 
tween formation and dissociation does lead to a very 
broad yield maximum in the vicinity of the expected 
range of initial temperatures at RHIC. This effect has 
been noted for the nominal expansion and initial volume 
parameters in Fig.|7| Here we see this effect persists for 
a variety of these parameters. 

In order to establish a lower limit for the final Be frac- 
tional abundance, we show by the dashed and dotted 
curves in Fig. || the effects of increasing the expansion 
rate, thus decreasing the lifetime of the QGP. We see that 
this effect is less than that generated by a change in initial 
volume, which controls the formation rate through the 
change in charm quark density. The different initial vol- 
umes are controlled by varying li. The set of 6 lines in the 
middle corresponds to li = 1 fm, (Vq = 52 fm'^), our stan- 
dard scenario, while the upper set of lines corresponds to 
the high density case li = 0.5 fm, {Vq = 26 fm"^), and 
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the lower set of lines represents the low density case, 
h = 2fm, {Vo = 105 fm^). 

In all of these calculations we have used an initial 
charm quark number Nc^ — 10, and for the range of 
temperatures expected at RHIC, this number will not 
change appreciably during the QGP lifetime. Eq. |l] then 
predicts that Nb^ will be exactly linear with the initial 
number of charm quarks, and our numerical results verify 
this. The calculated final Be yields utilizing the average 
initial A^co thus will correctly include the fluctuations ex- 
pected according to a binomial (or Poisson) distribution. 
This linear property is not evident in the results shown 
in Fig. 1^ in terms of the initial charm density, since there 
one is changing the density by changing the volume, and 
this also effects the lifetime of the QGP and through that 
the final Be abundance. 

Shown in Fig.|| is the fractional i?c-yield dependence 
on initial volume measured through variation of U for a 
range of initial temperatures. One sees that the temper- 
ature variation is less important than the initial volume 
effect, at least in the range of these parameters which we 
consider. For a fixed Tg, a change in the initial volume 
not only causes a inverse variation in the initial charm 
density pc oc but also influences the plasma life time 

1/2 

tl. Empirically it was found that tl oc is a good fit. 

— 1/2 

These two factors combine two give the roug h\yl^ ^'^ de- 
pendence of the final fractional Be meson yield shown in 
Fig.g. 
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FIG. 9. Be fractional yield at freeze out as function of the 
initial (width) parameter U [fm] for To = 0.5 GeV (dotted), 
0.4 GeV (dashed), and 0.3 GeV (solid). 



C. Sensitivity to Breakup and Formation Cross 
Sections 

The breakup cross section magnitude and shape is an 
essential part of our dynamical calculation, since in ad- 
dition to providing a time scale, it provides through de- 
tailed balance the relative magnitude of the formation 
cross section. The form in Eq. ^ has been used to es- 



timate the breakup rate of J/ip due to final state colli- 
sions with hadrons, via convolution with the gluon struc- 
ture function of the hadron ||]. Recently, there have 
been several additional attempts to model the hadron- 
quarkonium cross section JlGf , which has lead to results 
which typically are much larger and have a more rapid 
rise above threshold. Here we investigate the sensitivity 
of our prediction Be yields to such variations in the fun- 
damental cross sections. We show in Fig. ^ the change 
in the predicted Be yields at RHIC which follow if our 
breakup cross section in Eq. |^ is increased by a factor 
of 2. Also shown are corresponding results if the cross 
section is assumed to immediately rise to its maximum 
value just above threshold, with the same overall magni- 
tudes as before. All other model parameters are kept at 
the nominal values. One sees that in all cases the final Be 
yields are increased. One can understand this behavior as 
a combination of two effects. First, the detailed balance 
relation provides the relative magnitudes of formation 
and breakup rates. Second, an increase in the magnitude 
of the cross section just decreases the corresponding time 
scales, allowing the favored formation reaction to pro- 
ceed further toward completion. Thus the cross section 
we have utilized provides a conservative lower bound for 
our Be production estimates at RHIC via this new mech- 
anism. 



0.15 




0.05 



0.4 

Tg [GeV] 

FIG. 10. Be fractional yield at freeze-out for various the- 
oretical breakup cross sections, nominal (solid), 2x nominal 
(dotted), constant (dashed), 2x constant (dot-dashed). 



D. Relative excited state Be yield 

We have also calculated the ratio of 2S to IS-state 
Be yields within our model scenario. This is prompted 
by the observation that the corresponding ratio in the 
charmonium system at SPS may serve as a thermometer 
of the QGP phase 

As a first estimate of this ratio we use as a 2S dis- 
sociation cross section Eq.^ with e2s — 0.25 GeV. This 
of course is only a rough guess, since one should change 
the parameters to those corresponding to a 2S state, but 
this state has binding and energy levels which are very 
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marginal in terms of the constraints used for the vahdity 
of the cross section formula. The individual population 
equations can be solved independently, since both final 
bound state fractions are small enough that the source 
of b-quarks is not significantly decreased. 

Fig.0 shows the ratio of yields Bc(2S)/Bc(lS) at 
hadronization as a function of the initial temperature for 
the three different initial volumes we consider: h = 2 fm 
(dotted line), Ifm (dashed line), and 0.5 fm (solid line). 
It appears that this ratio is somewhat more sensitive to 
the initial temperature than the individual yields, vary- 
ing by at least a factor of two within the expected range 
for RHIC. 




0.2 0.4 0.6 0.8 1 

[GeV] 

FIG. 11. Ratio of Bc(2S) to -Bc(lS) abundance at 
freeze-out vs. initial temperature [GeV] for k — 2fm (dot- 
ted line), Ifm (dashed line), and 0.5 fm (solid line). 



We also find that this yield ratio is insensitive to ini- 
tial charm abundance and production, as charm density 
enters linearly in both the Bc{2S) and -Bc(lS) population 
equations. Also, the initial volume as shown in Fig.[ll| 
does not alter the abundance ratio significantly. Thus 
this yield ratio may allow one to draw conclusions about 
initial temperatures present in the QGP phase at RHIC, 
independent of the other parameters. To estimate the 
systematic uncertainty that would arise in such a proce- 
dure, we consider the effects of changing the cross sec- 
tions (formation and by detailed balance also breakup) 
in the following cases: (a) increasing the (2S) cross sec- 
tions by a factor of two, (b) increasing both 2S and Is 
cross sections by factor of two, and (c) decreasing the 
Bc{lS) and 5^(28) binding energies by 100 MeV. In all 
three cases the overall shape of the results as shown in 
Fig. |ll] remain nearly unchanged. However, we note an 
overall increase in the relative yield by a factor of two 
in case (a), practically no change in the result in case 
(b) and a decrease by a factor of 1/3 in case (c). The 
ratio Bc{2S) / Bc{lS) is thus primarily sensitive to initial 
temperature, but to be able to draw firm conclusions we 
need accurate relative 2S to IS cross sections, and a good 
understanding of in-plasma Be binding energy. 



V. SUMMARY AND CONCLUSION 

We have shown that this new mechanism of quarko- 
nium production in a deconfined medium predicts the 
minimum final state abundance of Be mesons at RHIC 
to be of the order of 5% per initial bb pair. This result is 
relatively insensitive to the initial temperature and vol- 
ume of the QGP, as well as to changes in the transverse 
expansion dynamics. There is a linear dependence on 
the charm quark abundance in the initial state, which in 
fact is the primary controlling factor in the final state Be 
yield. 

Fractional Be yields at the level of 5 x 10~^ signifi- 
cantly exceed expectations based on initial coherent one 
step production in individual nucleon-nucleon interac- 
tions, where a relative yield in range 10~* — 10~^ is 
expected. Such a small yield would not be observable 
at RHIC. Even the 500-1000 times greater multistep 
yield we obtain may pose considerable experimental chal- 
lenges. However, should such an experiment succeed to 
even roughly confirm these predictions for the fractional 
Be yield, we would have a convincing evidence for the 
mobility of charmed quarks over an extended space-time 
region in the dense phase. 

While in principle one could argue that incoherent Be 
formation could also occur in the hadronic phase in colli- 
sion of D-mesons with B-mesons, such a process requires 
localization in phase space of both these hadrons, which 
is highly unlikely. A calculation has been performed for 
the analogous case in the J/tp system, using D- meson 
interactions |l^. It was found that this mechanism is 
negligible even at LHC energies, except for possibly some 
observable effects in the tp' yield. 

We emphasize here that an essential element of this 
calculation relies on the assumption that colored heavy 
quarks will be subject to large energy loss processes in a 
plasma . This is necessary if these heavy quarks are 
to exist in a common region of the phase space volume. 
Without this stopping, heavy quarks are highly unlikely 
to remain for the required period of time inside the ther- 
mal deconfined phase. The details of this scenario cannot 
be fully justified at present, and considerable effort has 
to be vested to better understand the mechanisms which 
determine the initial phase space distribution of heavy 
quarks formed in nuclear collisions. 

Another issue is the final state reduction of Be mesons 
due to collisions with comoving hadrons. However, we 
would expect the smaller size of these bound states rela- 
tive to those extensively studied in the J/tp system will 
produce a much less significant effect in the Be system. 

Our work relies heavily on the presence of mobile 
charmed quarks: investigation of other bound state for- 
mation, such as J/tp, is currently underway. Initial re- 
sults indicate that while a significant fraction of charm 
will be in fact bound (at the 10% level) this is too little 
to significantly alter the results we have presented for Be 
yields. In fact the overall initial charm production uncer- 
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tainty is at least as large as this effect. One such source 
is the possible shadowing of gluons in nuclei, which could 
reduce the primary yield by amounts in the tens of per- 
cents Hi. 

If the experimental techniques enable observation of 
these predicted Be yields with significant efficiency, one 
has the possibility to embark upon a general study of 
both the vacuum and plasma properties of the Be sys- 
tem. Even should this prove not to be possible in the 
RHIC energy range, it will most certainly be easier at 
the Large Hadron Collider (LHC), where the collision 
energy is 30 times greater. Several bottom quark pairs 
will be produced in each central nuclear collision, along 
with literally hundreds of charm quark pairs, leading to 
expectations involving copious production of Be mesons. 
A study of relevant parameters for this production mech- 
anism at LHC is underway. 
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